A Correlation Between Surface Densities of Young Stellar Objects 
and Gas in Eight Nearby Molecular Clouds 
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ABSTRACT 

We report the discovery and characterization of a power law correlation between the local 
surface densities of Spitzer-identified, dusty young stellar objects and the column density of gas 
(as traced by near-IR extinction) in eight molecular clouds within 1 kpc and with 100 or more 
known YSOs. This correlation, which appears in data smoothed over size scales of ~1 pc, varies 
in quality from cloud to cloud; those clouds with tight correlations, MonR2 and Ophiuchus, 
are fit with power laws of slope 2.67 and 1.87, respectively. The spread in the correlation is 
attributed primarily to local gas disruption by stars that formed there or to the presence of very 
young sub-regions at the onset of star formation. We explore the ratio of the number of Class II 
to Class I sources, a proxy for the star formation age of a region, as a function of gas column 
density; this analysis reveals a declining Class II to Class I ratio with increasing column density. 
We show that the observed star-gas correlation is consistent with a star formation law where 
the star formation rate per area varies with the gas column density squared. We also propose 
a simple picture of thermal fragmentation of dense gas in an isothermal, self-gravitating layer 
as an explanation for the power law. Finally, we briefly compare the star gas correlation and 
its implied star formation law with other recent proposed of star formation laws at similar and 
larger size scales from nearby star forming regions. 

Subject headings: dust, extinction - infrared: stars - stars: formation 



1. Introduction 

It has been long recognized that the density of 
young stars varies by orders of magnitude be tween 
different star forming regions. Irlerbia ( 1962 ) com- 
pared the low density of young stars in the Tau- 
rus star forming region to the high density of 
young stars in the Orion nebula and speculated 



that these may form an evolutionary sequence, 
with the low density regions collapsing into denser 
structures. More recently, grou nd-based near-IR 
surveys of molecular clouds (e.g . Lada et all 1991: 
Strom. Strom, fc Merri ll 1993; Carpenter 2000; 



Porras et all 120031 lLada fc Lada I I2003J) identi- 
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tied embedded dense young stellar clusters and 
groups in molecular clouds as well as a more dif- 
fuse population of low mass stars, showing that 
molecular clouds contain both a dense "clustered" 
and a diffuse "distributed" population. Although 
clusters could easily be identified in these sur- 
veys by searching for peaks in the surface den- 
sity of sources, measurements of the number of 
distributed young stars in the clouds were domi- 
nated by uncertainties in the amount of contami- 
nation by field stars and the tre atment of variabl e 
extinction through each cloud (jCarpenterl |2000|) . 
This led to a number of questions regarding the 



number and ori gin of young low mass stars out- 
side of clusters (|Lada et al.l Il991[ iLi et al.l Il997 ; 
Carpenterll2000l ). What fraction of the cloud pop- 
ulation is found outside of clusters? Are they 
born in clusters and dynamically ejecte d to their 
observed positions (jEvans et al.l [2009)? Alter- 
natively, if they are formed near their current 
positions, why do molecular clouds exhibit large 
variations in the density of young stars? Are the 
variations in density related to variations in the 
Jeans length within clo uds ( Teixeira et al.l l2006t 



Gutermuth et al.l 120091 )? Better observational 



characterization of the relationship between the 
spatial distribution of YSOs and their natal cloud 
material is needed to understand the nature of the 
distributed component, constrain models of star 
formation in molecular clouds, and ascertain the 
underlying physics that det ermine the density and 
rate of star formation (e.g. Bate k Bonnellll2005c 



Krumholz k Tan 1 120071; iMvers 2009a) 



Sensitive Spitzer mid-IR and 24 /im surveys 
of nearby young clusters embedded in molecu- 
lar clouds have enabled complete characterization 
of the full range of stellar densities in a range 
of star foming environ ments (e.g. lAllen et al.l 
120071; lEvanset al1l2009l) . Young stellar objects 
(YSOs) in the clouds are identified by the pres- 
ence of IR-excess emission from circumstellar disks 
and/or envelopes. Several groups have devel- 
oped methods to use Spiteer-determined YSO 
spectral energy distributions to characterize the 
sources as protostars (Class 0, I, or flat-spectrum) 
or pre-main sequence st ars with disks (Class II 



or transition disk s ) (e.g. Gutermuth et al.l 120091 : 



Evans et al.l 120091 iRebull et al.l l2010h . By em- 



ploying strategies to minimize the contamina- 
tion by bac kground galaxies wi t h similar infrared 
color s (e.g. lHarvev et al.l 20061 ; iGutermuth et al 



20081 l2009t lEvans et al.l l2009l ). the Spit 



tzer sur- 
veys have yielded reliable, high spatial dynamic 
range distributions of YSOs in their natal clouds. 
These observations show that the clusters are 
peaks of a more extended population of young 
stars which tend to follow the morphol ogy of the 
filamentary molecular clouds ( j Allen et al.l [2007; 



Gutermuth et al.l 120091; lEvans et al.ll2009h . Thus, 



low density star forming regions are not limited to 
Taurus-like regions, but are also found in molec- 
ular c louds forming lar ge clusters and high mass 
stars (|Allen et al.ll2007l ). 



The wide range in observed YSO surface den- 
sities provides an opportunity to study how the 
density of young stars varies with the observed 
star formation efficiency and the properties of the 
co-spatial molecular gas. Recent work indicates 
that the star format ion efficiency increa ses with 
the stellar density. Evans et al.l (|2009l ) showed 
that higher surface density YSO clusterings tend 
to exhibit higher star formation efficiency (30%) 
than lower density surroundings (3-6%). In a 
Spitzer survey of YSOs in the W5 molecular cloud, 
where a recent generation of star formation is pro- 
ceeding on the peripheries o f the evacuated cavi - 



ties of adjacent HII regions, iKoenig et al.l (|2008l ) 



found star-formation efficiencies of >10-17% for 
high surface density clusterings and 3% for lower 
density distributed regions, respectively. The frac- 
tion of Spitzer identified YSOs outside of clusters 
depends strongly on the adopted functional defini- 
tion of stellar cluster, as there is no strong distinc- 
tion in source spacings evident in the nearby YSO 
spacing data to suggest distinct " clustered" and 
"dist ributed" star-forming modes (JBressert et al. 
2010)). The continuous nature of the overall YSO 



surface density distribution in nearby clouds mo- 
tivates a different approach whereby the local sur- 
face density of YSOs and the local column den- 
sity of gas are compared. This approach does not 
require a functional definition of clustered or dis- 
tributed star formation. 

A similar approach has been us ed with success 
i n obs ervations of external galaxies. iKennicutt et al 
(|l998l ) outlines previously published results con- 
cerning the validity of the Schmidt Law of star 
formation in external galaxies, namely Tisfr = 
AYi^„ a , where the Es refer to the observed surface 
densiti es of star form ation rate and gas respec- 
tively (|Schmidtlll959l ). For a limited range of den- 
sities and galaxy typ es these studies f ound N to 



be between 1 and 2. IKennicutt et al.l (|l998l ). us 



ing Ha emission as the star formation tracer, and 
CO and HI as the gas tracers for spiral galaxies, 
and Br7 emission and FIR dust emission as star 
and gas tracers respectively for infrared-selected 
starburst galaxies, show N = 1.45 ± 0.15. Note 
that these observations measure the integrated 
star formation (as traced by recently formed OB 
stars) across multiple molecular clouds, hence the 
extragalactic laws should not be expected to ap- 
ply directly to studies of low mass stars within 



individ ual molecular clou ds in our galaxy. That 
stated, Evans et al.l (J2009I ) compare T,sfr against 
T, gas for the seven nearby clouds surveyed by the 
Spitzer c2d Legacy survey against the Kennicutt- 
Schmidt Law, and find ^sfr to be ~20 times 
larger than predicted by the law for the measured 

^gas- 

In this contribution, we perform a direct com- 
parison between locally determined mass surface 
densities of YSOs and gas. This technique reveals 
a power law correlation between the YSO surface 
density and the gas surface density in all clouds, 
which we refer to hereafter as the star-gas corre- 
lation. Furthermore, using the ratio of protostars 
to pre-main sequence stars as an evolutionary in- 
dicatory, deviations from that correlation can be 
interpreted as regions of more or less time evo- 
lution in the star-gas correlation. We interpret 
these by adopting a star formation law of power 
law form and calculating the integral of the star 
formation rate over time vs. the remaining gas 
density. We present a simple analysis which shows 
that the higher density regions of a cloud produce 
stars which consume their local gas more rapidly, 
and the relation between log £* vs. log T, gas in- 
creasingly deviates from a power law with time, 
particularly at high gas densities. Nevertheless, 
the observed correlation can be reasonably well 
explained by an underlying Y*sfr — AY? gas star 
formation law. We examine the observed devi- 
ations from the correlation and argue that they 
result from a combination of gas dispersion in rich 
clusters and non-coevality within the molecular 
clouds. 

In Section [21 we summarize the Spitzer YSO 
lists compiled from prior observations, as well as 
the near-IR extinction maps (as a tracer for the 
molecular gas surface density). In Section [3l we 
demonstrate and characterize the observed corre- 
lations. In Section |4j we present some analysis 
of the correlation and its implications including: 
a simple evolutionary model to demonstrate the 
robustness of an underlying star formation law 
of power law form under the effects of early gas 
depletion; a simple thermal fragmentation model 
that predicts the central values of the ranges of 
power law index and normalizations observed; a 
brief comparison to the Schmidt-Kennicutt star 
formation law for galaxies. Finally, we summarize 
our findings in Section O 



2. Data and Methodology 

Numerous surveys of nearby star formation 
regions have been performed using Spitzer (e.g 
Gutermuth et al.l2003lMegeath et al. 2OO~il;|j0rgensen et al 



20061 lHarvev et al.l 120071: lAllen et al.ll2007t) . One 

of the most useful products of these observations 
is a largely flux-limited census of high probabil- 
ity young stellar object candidates as indicated 
by their excess mid-infrared emission relative to a 
photospheric spectral energy distribution. Spitzer 
data cannot be used to uniquely distinguish disk- 
less pre-main-sequence stars from field stars, but 
this limitation to the total YSO census is rela- 
tively minor for young, embedded star forming 
regions where few members have had sufficient 
time for thei r disks to disperse (10- 50% for ages 
<1 to 3 Mvr: iHernandez et al.l[2008l ). 

Table Q] lists the cloud surveys we have used 
for this analysis, the references for these surveys, 
and some general statistics for each. We selected 
large clouds with 100 YSO members or more, 
with a known range of YSO surface densities that 
is greater than one order of magnitude, and lo- 
cated within 1 kpc of the Sun, as these typically 
are complete for infrared-excess so urces down to 
rough ly 0.2 M Q at an age of 1 Myr ( Baraffe et alj 
19981 ). and to even lower masses (< 0.08 M Q ) 
for the closer clouds. However, the presence of 
bright sources as well as bright, structured nebu- 
losity near the centers of some clusters does sig- 
nificantly reduce the completeness relative to re- 
gions of lower YSO density. The most significant 
example of t his issue is the Trapeziu m area of 
Orion, where Meg eath et al.l (|in prep. I ) have per- 
formed a detailed comparison between the Spitzer 
and Chandra source lists to examine the issue of in- 
completeness. The completeness decreases rapidly 
in the inner regions of the Orion nebula: the com- 
parison has revealed that the Spitzer survey may 
identify 50% of the dusty YSOs averaged over the 
inner 0.5 parsec radius, and only identify 30% of 
the YSOs in the inner 0.2 pc. Other expected 
trouble spots include the centers of the MonR2 
and GGD 12-15 clusters in the MonR2 cloud, and 
the center of the Cep A cluster in the CepOB3 
cloud. In these regions of generally high YSO sur- 
face density, we are likely underestimating the true 
surface density of YSOs. 

The adopted datasets are relatively uniform 



in their observing parameters, but heteroge- 
neous in their data treatment, including different 
software pipelines used for mosaic construction, 
source extraction algorithms, photometric mea- 
surement techniques, and finally source classifi- 
cation schemes. While the results from the first 
three items are expected to be similar, it is not 
intuitive to ascertain the degree of equivalence of 
the YSO type classification schemes. Some inves- 
tigations into this issue have suggested that the 
relative numbers of Class I and II YSOs are found 



to va r y at the 10% level by technique (I Allen et ah 
20071 : iRebull et al.ll201Ct IGutermuth et al.ll2009t ): 



method dependent variation at that level should 
not strongly affect the results presented here. 

2.1. Large-scale Near-IR Extinction Map- 
ping 

To characterize the spatial distribution of gas 
column density in the target molecular clouds, we 
have adopted dust column density maps from the 
near-IR reddening of stars in the field of view, as- 
suming that the dust and gas column densities 
are proportional. The maps ar e constructed fol- 
lowing the method described by IGutermuth et al. 



(|2009h for the Orion A and B, MonR2, CepOB3, 
S140, and North America Nebula (NANeb here- 
after) clouds. For the c2d survey clouds presented 
(Ophiuchus, Perseus, and Serpens), the extinction 
maps delivered as part of that survey's enhanced 
data products have been adopted, after a compar- 
ison between ove rlapping regions o f the c 2d maps 
and the maps of IGutermuth et al.l (J2009I ) showed 
consistent results (see Appendix lAf . The former 
set of maps were generated at the Nyquist sam- 
pling resolution for the median n — 20 nearest 
neighbor distance in each region. The sampling 
resolutions of each map are tabulated in Table Q] 
though it is worth noting that the effective resolu- 
tion of these maps decays somewhat at higher col- 
umn density. The mean and maximum values after 
baseline subtraction, and the value subtracted as 
a baseline, are also inclu ded in the table. The em- 



pirical extinction law of iRieke fc Lebofskvl (|1985h 
was adopted to convert the H — K$ color excesses 
to A v (A V /E{H - K s } = 15.9). The gas column 
density maps are computed from the extinction 
maps by assuming a constant gas to dust ratio of 
A V /N{H 2 ) = 10~ 21 magnitudes-cm 2 . This ratio 
has been observed to vary by up to a factor of two 



(e.g. lWinston et al.ll2010h . however, such variation 
will have a modest effect on the star-gas correla- 
tion which is observed over more than an order of 
magnitude range in the gas column density. 

2.2. YSO Surface Density Estimation 

We have performed a n th nearest neighbor sur- 
face density analysis of the Spzteer-identified YSO 
distribution in each cloud, sim i lar to the analysis 
employed in IGutermuth et al.l (J2009I) . For these 
values, n was set to 11 - this value results in uncer- 
tainties of 33% in the s urface density estimation 
(|Casertano &: Hutlll985l ). We first obtained these 
nearest neighbors centered on each of the iden- 
tified YSOs; for a sample of q YSOs this would 
result in q surface density measurements centered 
on the YSOs. We also obtained nearest neigh- 
bor distances for a uniform rectangular grid of 
points in RA and Dec; this is the same procedure 
used to generate nearest n eighbor density maps in 



several recent papers (e.g. Gutcrmut h et al 



Roman-Zuniaa et al.1 120 08: Cha varria et al.l 
Jpreensen et al.l 120081 : IGutermuth et al.l 



2005; 



2008; 



2008; 



Evans et al.l 12009: G utermuth et al.l l2009h . To 



facilitate comparisons between the stellar density 
and the gas column density, the grid vertices we 
have adopted here are identical to the extinction 
map's pixel positions for each region. The com- 
bined analysis yields two surveys of the gas and 
stellar surface densities, one that follows the YSO 
spatial distribution and the other that is unbiased 
and spatially regular. In both cases, the resulting 
densities of YSOs per solid angle in the sky are 
converted to mass per square parsec using the dis- 
tances in Table [T] and adopting an average stellar 
mass of 0.5 M Q . 

3. An Observed Correlation in the Surface 
Density of YSOs and Dust 

Visual comparisons between the distribution of 
.Spiteer-identified YSOs relative to maps of gas 
structure in their natal molecular clouds have 
shown that there is considerable similarity be- 



tween them (e.g. iMeeeath et al.ll2004J: lAlleri et al 
20071 lEvans et al.l 120091: IGutermuth et al.l 120091) 



In Figure [TJ we show one of the best examples of 
cloud- YSO correlation, the MonR2 cloud. The red 
dots are the positions of Spiteer-identified YSOs, 
the grayscale image is the 2MASS-derived extinc- 



tion map, and the green line outlines the Spitzer 
mid-IR survey coverage. The YSOs appear to 
congregate in the regions of detectable extinction, 
peaking in density at or near the positions of high- 
est extinction. Figure [5] is the same figure, but for 
the CepOB3 cloud. Here, the apparent correla- 
tion is similar if one considers only the regions of 
significant extinction. However, there is an ex- 
posed young cluster, CepOB3b, in the northwest 
corner of the coverage where there are many YSOs 
present and very little local extinction; the gas in 
this region has likely been dispersed by winds and 
radiation from associated high mass stars. (See 
Figures |3][8] for similar figures for the rest of the 
clouds.) 

To quantify the apparent correlation, we plot 
(on a log-log scale) the YSO mass surface den- 
sity £* vs. the molecular gas mass column den- 
sity T, gas for MonR2, CepOB3, and six other star- 
forming clouds in Figures |H] & QUI utilizing the 
data and conversion factors described above. In 
Fig. [9] for YSO-centered sampling, red points rep- 
resent central YSOs classified as protostars (here- 
after Class I, though Class and flat spectrum 
sources are included), and grey points represent 
those that are classified as pre-main sequence stars 
with disks (hereafter Class II, though transition 
disks are included). The black hatched regions 
represent regions in the plot space where the ex- 
tinction is too low to be reliably measured in the 
near-IR at these resolutions (specifically Ay < 1) 
and where the surface density of stars is consistent 
with residual contamination from AGN (7 deg~ 2 ; 
Gutermuth et al.l I2009H H The bias toward sam- 



pling only in positions with YSOs is addressed in 
Fig. 1101 Here, plot points are sampled uniformly 
by position on the grid defined by the extinction 
map sampling, yielding an area-sampled, rather 
than YSO-sampled, view of the correlation. 

Regardless of the view used, every cloud exam- 
ined exhibits a positive power law trend in the 
surface density of stars and column density of gas. 
The MonR2 and Ophiuchus clouds stand out as 
having particularly cohesive loci relating these sur- 
face densities (Pearson coefficients of the log-log 
data are 0.87 and 0.83, respectively), thus we use 



1 This value is highe r than the value reported for c2d (e.g. 
lHarvev et al,l 120071 ). but we apply it across all clouds for 
consistency in treatment. 



their correlation data to characterize the trend in 
detail. Using the sa me fitting method used by 
iFlahertv et~atl (J2007FF1 . we fit lines to the log £* vs 
log E gQS data for MonR2 and Ophiuchus, finding 
them well fit with power law indexes of 2.67 ±0.02 
and 1.87 ± 0.03, respectively. Uncertainties in 
these values are simply the intrinsic slope uncer- 
tainty from the fit; we estimate a systematic un- 
certainty of < 0.04 based on a comparative anal- 
ysis of Ophiuchus using both our extinction map- 
ping technique, and that of c2d (see Appendix lAl. 
The fits are overlaid as polygons (green/steeper 
and blue/shallower for MonR2 and Ophiuchus, re- 
spectively) on the measurements for all clouds in 
Figs.[H]&[TU]to demonstrate that the correlations 
for the other clouds are largely consistent with a 
power law of index ^2. However, most clouds have 
some larger spread in the observed correlation rel- 
ative to MonR2 and Ophiuchus (~1 dex spread), 
and this fact deserves closer examination. We have 
included fit results for all of the clouds in Table [51 
for reference. 

The more complex relationships seen in the 
other clouds appear related to evolutionary dif- 
ferences in the star-gas correlation among differ- 
ent structures within these clouds. For exam- 
ple, the main cluster core of the IC 348 region 
is well known to be a more evolv ed population 



of YSOs, with evolved disks f e.g. iMuench et al.. 
2007), a low disk fraction (e.g. iLada et alll2006l) . 



a hi gh ratio (27) of Class I I to Class I sources 



20091) . and low mean ex- 
tinction relative to other young clusters (e.g. 



(e.g. iGutermuth et al. 



J0rgensen et all2008t [Gutermuth et~al~ll2009l ). All 
of these observations indicate an evolved YSO 
population and thus significant time for the YSOs 
and the intermediate mass central members to 
have dispersed the local gas via feedback pro- 
cesses. In the YSO vs. gas surface densities di- 
agram for Perseus, the subset of points that are 
significantly over-dense in stars relative to the 
gas density predominantly belong to the IC 348 
exposed cluster core. Perseus also has several 
regions that are under-dense in stars relative to 
gas. These regions are typically located in the 



2 The lFlahertv et al. (2007) algorithm is an IDL implementa- 
tion of x 2 minimization of the two dimensional perpendic- 
ular offsets from a line fit, thereby implicitly incorporating 
uncertainty in both axis quantitites, following the analysis 
of iPress et all (|1992T ). 



well-studied dark clouds in Perseus, such as B5 
and L1448, where the YSOs in these objects are 
predo minantly confirmed Cla s s and Clas s I ob- 
jects (|j0rgensen et all 120081: lEnoch et all l2008h 
with very few more evolved Class II sources. This 
indicates relative youth when compared to more 
typical embedded clusters like NGC 1333, which 
has a more standard ratio of Class II to Class I 
sources (3.7: iGutermuth et al.ll2009l ). Ultimately, 
this study reveals that the Perseus cloud is an 
outlier relative to other nearby molecular clouds 
because of its heterogeneous composition of evo- 
lutionary states among its star-forming sites. 

The Cep OB3 giant molecular cloud (GMC) 
appears to be a simpler (yet more numerous 
in YSOs) analogue of Perseus in that it has a 
subset of YSOs that deviate from the observed 
trend, namely the large, exposed Cep OB3b clus- 
ter consisting of two sub -clusters, Cep O B3b- 



east and Cep OB3b-west (jAllen et al.l lin prepl ) 



Cep OB3b-east has an abnorma lly large cur- 
rent star-formation efficiency f80%: lGetman et al. 



2009). The cluster is found to be adjacent to two 
massive molecular clumps indicatin g that most of 



the natal gas has been d ispersed (jGetman et al 



2009c lAllen et alJIin prep.l ) . Interestingly, the sim- 
ilarity of heterogeneous evolutionary states in Ce- 
pOB3 and Perseus may continue into the very 
young regime as well; there are some regions in the 
Cep OB3 molecular cloud which are overdense in 
gas relative to YSOs hinted at in the area-sampled 
analysis shown in Fig. [TO] 

Given these qualitative regional evolutionary 
differences, we investigate the degree that YSO 
evolution may trace differences in the S* vs. S gas 
analysis. We visualize the differing Class II to 
Class I ratios (Nn/Nr) of regions above or be- 
low the power law trend by assuming that the 
power law index is always 2 and computing the im- 
plied offset coefficient for the power law for each 
YSO individually. These values, £*/£ 2 QS , then 
may be composed into histograms separated by 
YSO evolutionary class. We have produced such 
histograms for all of the clouds considered here 
(see Figure [TT]) . The black histograms are the 
Class II YSO values, and the red histograms are 
the same, but for the Class I YSOs scaled up by 
3.7, the median rat io of Class II to Class I sources 
in nearby clusters (JGutermuth et al.l 120091 ) . Both 
have been normalized to the peak value among 



both histograms. The green histogram plots the 
Class II YSO data, but now scaled to the Class I 
data's normalization and divided by 30, as < 
1/30 Class II objects may be misidentified as 
Class I YSOs due to an edge-on disk orientation 
(JGutermuth et al.l 120091 ) . Internally coeval clouds 
like MonR2 and Ophiuchus appear to have very 
similar distributions for both YSO classes. In con- 
trast, those that are more heterogeneous in their 
evolution have a clear systematic trend whereby 
there is an excess in the Class I distribution where 
the calibration value is low (gas-rich or YSO-poor) 
and an excess in the Class II distribution where the 
calibration value is high (gas-poor or YSO-rich). 

Empirically, we define Y^/Y? gaa = 3 x 10~ 4 
and 5 x 10~ 3 pc 2 Mg 1 as boundaries that bracket 
the distributions of those clouds with clear cor- 
relations and set apart the extrema of local evo- 
lution that appear commonly in the more het- 
erogeneous clouds. Note that we have reduced 
over three orders of magnitude variation in £* to 
an approximately one order of magnitude varia- 
tion in S»/S 2 QS . S*/Sg as < 3 x 10~ 4 indicates 
a relative overdensity in gas, and the stars that 
are present there are predominatly of the younger 
type, Class I. We infer that the the onset of star 
formation was delayed in these regions, resulting 
in the presence of few YSOs and a large reser- 
voir of molecular gas. On the other extreme, 
£»/£ 2 as > 5 x 10 -3 suggests stellar overdensity 
relative to gas, and the stars found in this regime 
are nearly exclusively Class II YSOs. Thus we in- 
fer that these are regions where star-forming gas 
has been largely depleted by some combination of 
conversion into stellar mass or ejection of gas via 
radiative processes, such as photoevaporation, or 
winds from the recently formed stars. We sum- 
mariuze the median X^/Xl 2 ,^ values by YSO class 
for each cloud in Table [3] 

A majority of the YSO membership in each 
cloud (54%-85%) falls in the range of £»/ S gas 
between 3 x 10~ 4 and 5 x 10~ 3 pc 2 Mq 1 . This 
excludes regions of each cloud that have mostly 
dispersed their gas or appear to be just at the on- 
set of star formation. Such a selection has been 
applied to make Figure Q21 the Class II to Class I 
ratio of sources as a function of gas column den- 
sity for each of the clouds. This analysis yields a 
diagnostic of the evolutionary state of the molecu- 
lar clouds at different gas column densities. Even 



within the restricted range of £*/£g as , there are 
significant variations in the ratio of Class II to 
Class I objects. Overall, most clouds show a trend 
of decreasing Class II to Class I ratio with increas- 
ing gas column density. This may suggest that the 
higher column density regions are younger than 
those with low column density; however, there 
are several other possible interpretations for this 
trend. For example, the protostar phase may last 
longer in high density regions, or protostars in 
these regions may be less prone to oscillate be- 



tween Class II and Class I SEDs (jDunham et al 
120101 ). Alternatively, the regions of low column 
density may be contaminated by Class II objects 
that have migrated from clusters. Gas dispersal 
may result in the high fraction of Class II at low 
gas densities; this may be especially important for 
the Cep OB3b, S140 and Perseus clouds. Finally, 
we cannot rule out observational biases. Although 
our adopted classification criteria were designed to 
minimize the influence of extinction, high extinc- 
tion may result in the mis-classification of some 
Class II YSOs as Class I protostars. 

4. Discussion 

The Class II YSO evolutionary stage las ts be- 
tween 0.5 and 5 Myr ( Hernandez et al.ll20 08). and 



the majority of the IR-excess sources included in 
this analysis are of that class. Therefore, the fact 
that the above power law trend is observed at all 
suggests a close link between the YSO spatial dis- 
tribution and the gas structure on the timescale 
of ~2 Myr for the ~1 pc size scales probed here. 
While this is indeed a profound result, it is still 
in some sense unconstraining. The connection 
between the gas and the YSOs may in part re- 
sult from the observed kinematic connection be- 
tween dense cores and the surrounding molecu- 
lar gas. Specifically, millimeter-line observations 
show that the the velocity dispersion between 
the cores and the surrounding, moderate density 
gas is small relative to the velocity dispersion in- 
ferred from the observed linewidth and the ob- 
serve d velocity gradien t s in the bulk gas motions 
(e.s. IWalsh et all 12004 l2007t lAdams et al.l l2006t 
Kirk et alJl2007t JJeirgensen et al.ll2008l ). However , 
the YSOs and their progenitor gas could simply 
co-move in aggregate around the larger structure 
of the cloud, perhaps even dynamically interact- 
ing with other such aggregations, and still pre- 



serve the observed connection between the YSO 
and gas distribution. Thus, a dynamical picture 
for the cloud and the structures within the cloud 
is not ruled out by the star-gas correlation. 

Under the simplifying assumption that there is 
relatively little motion between the YSOs or the 
gas, we present a simple model which adopts a star 
formation law where the star formation rate per 
area follows a power law of the local instantaneous 
gas surface densities. This model demonstrates 
that the likely gas depletion that has occured 
to generate the observed star-gas correlation still 
leaves behind a power law of similar index to the 
primordial relationship. We then establish that 
the observed star-gas correlations and the star for- 
mation law have a power law index ~2. Next, we 
demonstrate that simple thermal fragmentation of 
an isothermal, self-gravitating, modulated layer of 
gas can produce a spatial distribution of gravita- 
tionally unstable fragments where the surface den- 
sity of fragments scales as the gas density squared. 
Finally, we briefly compare the reported star-gas 
correlation to similar recent analyses in the liter- 
ature. 

4.1. Comparing Star Formation Laws to 

the E* vs. Ilg as Relationship 

As discussed in Section [TJ several empirical star 
formation laws have been proposed in the liter- 
ature; these laws relate the inferred star forma- 
tion rate per area to the current column density 
of gas, and typically they are derived over large 
(>100 pc) size scale measurements. We examine 
here whether the observed ~1 pc scale £* vs. T, gas 
correlation reported in Section [3] is consistent with 
a simple star formation law where the star forma- 
tion rate per area varies linearly or as a power of 
the current gas column density. 

In order to consider the likelihood that the re- 
ported power law correlation in observed quanti- 
ties, ie. the integrated number of YSOs and the 
remaining amount of gas, relates at all to an un- 
derlying star formation law of power law form, we 
must consider whether the more efficient creation 
of stars at higher gas column density significantly 
affects the underlying gas distribution and there- 
fore results in an observed S* vs. T* gas distribution 
that deviates from the power law index of the star 
formation law itself. To this end, we present a sim- 
ple model of gas to stellar mass conversion. Gas 



mass is accreted onto stars at a rate governed by 
an assumed star formation law with an efficiency 
that accounts for the ejection of gas mass from the 
system via protostellar outflows. 

The nature of this model demands several un- 
derlying assumptions to both achieve adequate 
simplicity and offer a relatively robust examina- 
tion of how the observables behave with time. The 
first assumption is that the stars formed do not 
move significantly from their birth sites to within 
the typical resolution of our measurements (~1 pc 
scale). The observations that spatially proximate 
prestellar cores are likely formed with low relative 
velocities (|Walsh et al.ll2007t iMuench et aj]|2007t ) 
suggest that the typical velocity dispersion of of 
the YSOs that emerge from the cores is small 
(~0.4 km s^ 1 ). Furthermore, later in the star 
formation process, relatively uniform surface den- 
sity distribution s are observed in many embedded 
YSO groupings (JGutermuth et al.l 120091 ) . suggest- 
ing relatively little dynamical evolution. The sec- 
ond assumption is that there is no molecular gas 
flowing into or within the molecular cloud at the 
~1 pc scale; gas is only ejected from the system 
by the formation of stars, or converted into stellar 
mass. This assumption is adopted for simplicity; 
it is demonstrably false at scales <1 pc, where 
systematic gas infall motions have been observed 
in cluster- forming clumps (e.g. IWilliams fc Mversl 
f 999; Walsh ct al. 2006), and on larger scales, the 
global freefall time of a typical GMC is ~5 Myr, 
though turbulent support likely inhibits global col- 
lapse to some degree. Finally, we assume that 
there is no effect on the gas distribution by large 
scale effects from phenomena such as feedback 
from nearby high mass stars, supernovae, or galac- 
tic motions. These processes only affect the cloud 
on longer timescales (>10 Myr) than we are con- 
sidering, with the exception of high mass star feed- 
back that we ignore for simplicity. 

We begin by adopting a general star formation 
law where the star formation rate per area shows a 
power law dependence on the gas column density: 



dt 



ckljg as yx, y, t) 



(1) 



where £*(a;, y,t) is the mass surface density of 
YSOs as a function of position in the sky and 
time, k is a constant that is defined in equation [5] 
below, and S„ as (x, y, t) is the mass column den- 



sity of molecular gas as a function of position and 
time. We have assumed that the formation of a 
star depletes the cloud by M+/c where M* is the 
mass of the star created and c is the mass conver- 
sion efficiency which takes into account the mass of 
gas ejected from the cloud by feedback during star 
formation (M e j ecte d — M*/c — M*). Here we are 
assuming that the fraction of the mass in the out- 
flow ejected from the cloud is constant throughout 
the cloud, but in fact the fraction of mass that es- 
capes the cloud's gravitational potential may vary 
as a function of local gas column density and be 
affected further by the nearby gas configuration 
at >1 pc scales. Furthermore, it is worth noting 
that in this formulation we are not considering in- 
dividual stars, thus the entire analysis, including 
the value of c, is averaged over the stellar initial 
mass function (IMF), regardless of its form. A pre- 
stellar core to star mass conversion efficiency that 
is constant over the stellar mass range is consistent 
with recent characteriza tions of the mass function 
of pr e -stellar cores (e.g. lAlves. Lombardi fc Lada 



20071 : lAndreet al.l l2010h . They find a value for 
this efficiency of 0.3 ±0.1. However, this does not 
speak to the amount of unaccreted matter that is 
ejected from the cloud entirely; thus we consider 
0.3 as a lower limit on the value of c. 

If we assume that no additional gas is flowing 
into or within the cloud as the stars form, then 
the rate of depletion of the gas mass is given by: 



dT, gas (x,y,t) 
dt 



= -kZ gas (x,y,t) a (2) 



In the case that a^l, the solution for ~E gas (x, y, t), 
the gas density remaining in the cloud, is: 



Ega.(x,y,t) = £»„(&, W,0)(^- + If (3) 

to 

where, T, gas (x, y, 0) is the column density of gas at 
the onset of star formation, /3 is a constant, and 
to is the timescale for the gas to be depleted by a 
factor of 2^. Specifically: 



P 



I- a 



and 



fe(a-l) 



^>gas(x,y, 0) 1 



(4) 



(5) 



Note that to depends on the initial column density 
of gas for a > 1; in this case regions with higher 
column densities will be depleted more rapidly. 
The amount of mass converted into stars is then 
the mass of the depleted gas times the mass con- 
version efficiency, c: 



T,*(x,y,t) = c[Y gas ( x ,y,0) -T, gas (x,y,t)] (6) 

Substituting equation Eqn.[3Jinto Eqn[Sl we obtain 
the relationship: 



Z±(x,y,t) = cZ gas (x,y,0)[l - (- + 1)"]. (7) 

to 

In the case that t << to, we find that £*(#, y, t) 
has the same power law dependence as the star 
formation rate: 



£*(#,£/, *) = cY< gas (x,y, 0) a kt 



(8) 



In other words, the surface density of formed stars 
is proportional to the surface density of initial gas 
to the a power. In this limit, the observed gas 
density can be approximated as H gas (x,y,t) = 
^gas{x,y,0) and the following approximation for 
the star formation rate can be used: 



<9£*(x,2/,t) E*(a;,2/,t) 



dt 



t 



ckZ gas (x,y,t) a (9) 



In the case that the star formation rate is pro- 
portional to the gas density, i.e. a = 1 (even 
though our data cannot be reproduced by such a 
law, as we next show), we have a different solution 
to equation [2] for the gas remaining, H gas (x,y,t): 



At this point, we have argued that a power law 
observed in E* vs. E gas is indicative of a star for- 
mation law of similar functional form. Therefore, 
it is interesting to examine the a = 2 case, as sug- 
gested by the analysis from Section [3J 

Combining Eqn. [3] and [5] under the assumption 
that a — 2 allows us to find a useful expression for 
the quantity E*/E;L S computed in Section [31 



_£*_ 

Y 2 

gas 



(x,y,t) = ckt(l + -) (13) 

to 



Thus for a — 2, the histograms presented in 
Fig. [IT] can be thought of as histograms of nor- 
malized star counts for Class I and II YSOs 
vs. cktsF(x,y)(l + ^-(#(3/)) at the current 
mean age of the YSOs at a given position, 
t = tsF- Using this result and the observed evo- 
lutionary phase ratios of the YSO population, 
Nn/Ni vs. Y, ga s{x,y,t S F) (Fig. US}, we can con- 
strain the constant k and thus the gas depletion 
time to (E 9as ). Then we can directly compare our 
data to the gas depletion model. 

Here we measure tsF(x, y) by computing the ra- 
tio of the surface density of stars formed divided 
by the star formation rate surface density. We as- 
sume that the star formation rate surface density 
is constant in time and that it is approximated by 
the ratio of the surface density of Class I objects 
over the lifetime of t hat evolutionary p hase, E_r/tj, 
where t 7 ~0.5 Myr (|Evans et alJl2009h . Thus: 



tsF(x, y) = £*(£, y, t SF ) x ti/T, I {x, y) (14) 

?,*(x,y,tsF) = Bj(x,y) + Y II (x,y,t SF ) in this 
analysis, thus we can simplify to: 



-kt 



Y gas (x, y, t) = Y gas ( x , y, 0)e 
and for £*(x, y, t) we find the solution: 



£*(a;,2/,*) = cE oas (a;, y,0)[l - e 



-fc/i 



which becomes, in the limit of kt << 1: 



(10) 



(11) 



£*(x,y,t) = cY, gas (x,y, 0)kt (12) 

Eqn. [TU thus simplifies to Eqn. [5] for the a = 1 
case, under a similar assumption of small t. 



N n (x,y,t SF ) . , . 

t S F(x,y) = ( Nj ^ y) +l)xtj (15) 

As shown in section [3j we have computed 
Nu/Nj for each cloud as a function of gas column 
density, noting that some clouds have a down- 
ward trend, while others are fairly stable over the 
whole range of density values. The more extreme 
variations occur in clouds with wider spreads in 
£*/£;L s , thus we assume that those clouds have 
some contamination of lower gas column density 
regions by older YSOs (via their motion away 



from their birth sites) or that a large fraction of 
the primordial gas in those regions has been re- 
cently ejected by feedback from nearby, high mass 
stars. Given these complications, we have chosen 
to characterize the trend in Nu/Ni vs. S gQS from 
the relatively uniform MonR2 data, yielding: 



jj-(x,y,tsF) 



25 x T. gas (x,y,tsF) 



-0.4 



(16) 



The observed data are largely confined to the 
range 20 < £ 9as (x, y,t SF ) < 300 M Q pc~ 2 ; 
evaluating eqn. \TE\ at these fiducial points yields 
Nu/Ni = 7.5 and 3, respectively, and t$F — 4.3 
and 1.8 Myr, respectively. We now use these re- 
sults to set limits on the value of k, and thus to- 

Since tsF and to depend on quite different pow- 
ers of T, gas , the case of tsF << to likely occurs at 
low initial gas column densities (this will be veri- 
fied later), thus the higher value, tsF — 4.3 Myr, 
is applicable. Assuming c = 0.3 (the lower limit), 
then 2.3 x 10~ 4 < k < 3.9 x 10" 3 Myr' 1 pc 2 Mg 1 . 
Therefore, at the corresponding value of S gas 
(20 Mq pc~ 2 ), where the depletion time is long 
and thus 'E gas (x, y, t) = "E gas (x, y, 0), we find that 
the depletion time is 13 > t > 220 MyJE This 
range is 4 to 50 times larger than tsF = 4.3 Myr 
and thus consistent with the tsF << ^o assump- 
tion. Assuming that the spread in the MonR2 
£*/£ 2 s values are dominated by data uncer- 
taintjQ, we adopt the log-scaled mean value for 
the a = 2 case, k = 9.5 x 10~ 4 Myr^ 1 pc 2 Mq 1 . 

With k adopted, we now determine the gas col- 
umn density where to — tsF, as this is the regime 
where significant gas depletion is expected in the 
data and thus where deviations from a power 
law star-gas correlation would be detectable. For 
a = 2, to is the time at which £ gas (;r,y, to) = 
0.5 T, gas (x,y,0). Therefore: 



3 Note that to would increase (as k would decrease) linearly 
if c is larger, and c < 1 by its definition above. 

4 Our stellar surface densities are n = 11 nearest neighbor 
densities, which arc inherently uncertain at the 33% level 
(Cascrtano & Hut 1985)), and our gas surface densities are 
uncertain at 0.5-4 Ay, with higher uncertainty at higher 
surface densities. These uncertainty sources result in la 
uncertainties of ~0.3 dex in E t /S^ as values, thus they are 
sufficient to account for the ±0.6 dex range in those values 
in MonR2, assuming that they arc log- normally distributed 
and that our adopted limiting values isolate the central ±2cr 
range. 



1 

2k' 



to = 7^7Eg as (x,j/,to) 1 



(17) 



Setting to = tsF and thus this equation equal 
to Eqn. [15] and solving numerically, we find 



^gas(x,V,tSF = t a 



300M Q pc" 2 . Thus gas 



depletion only becomes noticeably relevant at gas 
column densities on the high end of the observed 
range. 

In summary, the model closely follows the ob- 
served power law correlation between star forma- 
tion rate surface density and gas column density. 
Only at the highest observed column densities 
does the model predict significant deviations from 
a power law: these may be undetected due to the 
scatter in the observed correlation. 

At this point, it is useful to visualize the time 
evolution of the simple gas depletion model and 
compare it with our observations. In Fig. 1131 
we plot the realizations of Equation [7] for sev- 
eral values of t in the range of 200 years to 
6.3 Myr and a range of S 90S (t = 0) from 20 to 
2000 Mq pc -2 , assuming a — 2, c = 0.5, and 
k = 9.5 x 10~ 4 Myr -1 pc 2 Mq 1 . These realiza- 
tions are isochrones showing the locus in the £* vs. 
£ gas diagram for an ensemble of co-eval molecular 
cloud parcels with different initial values of E ffas . 
We also plot orthogonal tracks that trace the evo- 
lution in £* vs. Egas for parcels of a molecular 
cloud as the gas of a given £ gas (t = 0) is converted 
into stars. The polygonal regions that demarcate 
the observed data limits for MonR2 (green) and 
Ophiuchus (blue) are overplotted for comparison. 
At early t, the model follows a power law of index 
~2. As we approach typical star formation region 
ages of ~1 Myr, the highest gas density evolu- 
tionary tracks curve toward decreasing gas den- 
sities, resulting in local steepening of the model 
isochrone. This is a direct result of the varying 
relative rate of gas depletion at differing column 
densities. For the lowest initial gas densities, the 
evolutionary tracks follow an almost vertical tra- 
jectory as stars are formed at very low efficiency. 
In contrast, the higher initial gas densities form 
stars very efficiently and thus the gas is quickly 
depleted relative to the star formation time, tsF- 
We can easily rewrite Eqn. [13] to describe the star 
formation efficiency at t = tsF as: 



10 



2.J*. 



-(x,y,tsF) 



^gas(x,y,tsF) 



N*=Nj 



,(x,y,t 



SF) 



l/{ckt SF (l 

(18) 



At the T, gas fiducial points of 20 and 300 M & pc 



adopted above, 



-(x,y,t S F) 



2.3% and 
to reported 



26% respectively. These are similar 
values for "di stributed" and " c lustered" star for- 
mation (e.g. ICarpenterl l2000t lAllen et al.l 12007 : 



mation (e.g . 
lEvans et al.l 



20091) . In th is analysis, as is claimed 
by iBressert et al.l (|2010l) . we see a smooth varia- 
tion from low to high column density star forma- 
tion. We add to that analysis that the most likely 
progenitor of the smooth variation in stellar densi- 
ties is the smooth variation of gas column densities 
in star forming clouds. 

In Fig. [T4J we investigate different values of a, 
as well as the possibility of a threshold in den- 
sity for star formation, and compare the results to 
the fits to the MonR2 and Ophiuchus clouds. We 
show Equation [7] models for a = 1 (top left panel) 
and a = 1.4 (top right panel), with no thresh- 
old imposed, and with appropriate values of k {a- 
dependent) to approximately fit the MonR2 and 
Ophiuchus data (green and blue polygons, respec- 
tively) : neither provides a good fit to the data, un- 
like that in Fig[T3J We also implement a threshold 
(bottom two panels) such that the star formation 
ceases when T, gas (x, y,t) is less than a value of 
100 Mq pc -2 . A less extreme threshold has been 
reported in a different analysis of the Spitzer c2d 
and Gould Belt Legacy survey data, where proto- 
stars alone are used to compute star fo rmation sur- 



face densities ( Heiderman et al.ll2010r ). Our anal- 



ysis is not consistent with a star formation thresh- 
old at any gas column density in the range that we 
have surveyed. 

4.2. A Simple Model of Jeans Fragmenta- 
tion 

For MonR2 and Ophiuchus, the power law in- 
dices in the relation between stellar and gas sur- 
face density are close to 2. This section shows 
that such an exponent is expected from a sim- 
ple model of Jeans fragmentation in an isothermal 
self-gravitating layer. 

If one initial Jeans mass makes one star, then 
the number surface density of stars in a layer 
equals the number surface density of Jeans masses, 



"IfThe Jeans mass in a self-gravitating layer is 



M, 



Ba 4 



i-'gas Ct 



(19) 



(20) 



where B is a constant ( 4.67 for an iso thermal, self- 
gravitating layer of gas; lLarsonlll985l ). a is the gas 



velocity dispersion, S gas is the gas mass surface 
density, and G is the gravitational constant. 
Then the number surface density of Jeans masses 
can also be written as 



Nj = 



or substituting from[20l 



Nj = -( 



S„ a 



gas 



Mi 



_2 > 



(21) 



(22) 



and substituting from ll91 multiplying both sides of 
the equation by a mean mass per star of 0.5 Mg|j, 
and converting to mass surface density units of 
Mq pc -2 , with an assumed a — 0.2 km s — 1 , the 
sound speed at 10 K: 



XL = 1.34 x 10 



-3y2 
gas 



(23) 



Thus the power law of index 2 is predicted 
by simple thermal fragmentation of an isothermal 
self-gravitating layer. Furthermore, the value of 
1.34 x 10~ 3 is in the middle (in log-space) of our 
observed range of £*/£ 2 as from section|3l and rep- 
resents a mean value for cfci(l + t$F /to) in the 
evolutionary model above. We note that there is 
no explicit time dependence in this model, thus 
although it gives the correct star-gas correlation, 
it does not predict a rate of star formation. 

More generally, one must consider cloud layers 
that are non-uniform, spanning a wide range of 
column densities in order to replicate the trends 



5 Note that the eventual mass of a given star is explicitly as- 
sumed to be uncorrelated with the local Jeans mass. This 
is consistent with some characterizations of the pre-stcllar 
core to stellar mass conversion process that include poten- 
tial accretion of gas from the large r-scale environmental 
distribution (e.g. lMversll2008l , l2009bt) . 
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we report here, which are found within each star- 
forming cloud individually. An example of a hori- 
zontally nonuniform planar system is the density - 
modulated isothermal layer (|Schmid-Burgklll967t) , 
or of a plana r system of isothermal cylinders 
( Ostrikerlll964 ) can be expected to give a similar 
relation between stellar and gas surface densities. 
Some of the clouds studied in this paper may have 
such a flattened geometry; their extinction maps 
indicate that "hub-filament" structure is expected 
if a clumpy medium is co mpressed into a modu- 
lated self-gravitating layer (|Mversll2009al ). 

If this picture of fragmentation in a layer ap- 
plies to some of the clouds considered here, the 
results of Section l4~Tl indicate that evolutionary ef- 
fects probably do not change the relation between 
star and gas surface density very much, over the 
likely star-forming life of the cloud. 



4.3. 



Empirical Comparison to Cloud-Scale 
Values, the Kennicutt- Schmidt Law, 
and Recent Parsec-scale Work 



As already noted, Evans et al.l (|2009f ) compared 
the overall star formation surface density vs. gas 
column density for each of the clouds in the c2d 
survey, as well as an average of all of them, for 
comparison with the Kennicutt- Schmidt law for 
galaxies. They find that the star formation rate 
density predicted by the law falls a factor of 20 
below their observed data. In Figure [15] we have 
overlaid the correlation reported here for even 
smalle r scales (0.3 t o 2.0 pc) on our reproduction 
of the Evans et al.l (120091) figure. Our correlation 
overlays the Evans et al.l ([20091 ) cloud-averaged 
data well at intermediate surface densities, and 
the mean value for all clouds falls quite close to 
our correlation. The difference is only clear where 
our trend extends to higher column densities. Be- 
cause we are able to probe a wider range in gas 
column density by taking measurements at smaller 
size scales, it is not surprising that such a discrep- 
ancy could be uncover ed. Regardle s s, our analy- 
sis agrees with that of Evans et al.l ([20091 ): stud- 
ies of star formation rate density versus gas den- 
sity for nearby clouds show much greater star for- 
mation rate for the local gas density than extra- 
galactic correlation s would suggest. As noted by 
Evans et al.l ( 20091 ) and others, the measurements 
used to derive most star formation law character- 
izations vary widely in the size scales and time 



scales probed, and these variations almost cer- 
tainly play a role in the apparent discrepancies. 
Indeed, recent work to characterize the role of 
measurement size scale in characterizing the cor- 
relation between star formation rate surface den- 
sities and gas surface densities in nearby galax- 
ies has revealed considerable change in the power 
law index with the measurement size scale used 



(|Liuet al.H2011l ). 

Other recent efforts have examined the star for- 
mation rate and gas surface density correlation in 
nearby star-formation regions via similar sorts of 
data (ie. YSO star counts and near-IR extinction 
maps) a nd some ov erlap in the specific clouds con- 
sidered. lLada et alj (|2010D estimated the star for- 
mation rate vs. cloud mass for a sample of 11 
nearby clouds. Within this cloud sample, they 
noted that there was a wide dispersion in the es- 
timated star formation rate per cloud mass. They 
then demonstrated that the observed dispersion 
within their cloud sample was minimized if they 
only considered the cloud mass at column densities 
above Ay — 7 (or 116 M Q pc -2 ). They proposed 
that this corresponded to a density threshold for 
star formation; for densities above this threshold 
the star formation rate varied linearly with gas 
mass, below the threshold the star formation rate 
was negligible. 

A similar threshold (129 ± 14 M Q pc -2 ) 
was invoked by the su b-cloud-scale analysis by 
Heiderman et al.l (I2010l). a follow-up to the cloud- 



p to the ( 
1 d2009h 



scale examinat ion by Evans e t alj . 
tioned above. Irleiderman et al.l ( 20101 ) used the 
numbers of Spiteer-identified Class I and so-called 
"fiat spectrum" YSOs (ie. they excluded the 
longer-lived Class II YSOs) within various ex- 
tinction contour intervals to probe the correla- 
tion in 20 nearby molecular clouds. Since their 
Spitzer data did not include clouds forming mas- 
sive stars, they did not have a direct measure- 
ment of YSO densities in high mass star forming 
regions. Thus they added star formation rates de- 
rived from IRAS far-IR luminosities and mean gas 
column densities derived from HCN line emission 
for a set of more distant and presumably more ac- 
tive star-forming regions. We have overlaid their 
final bent power law best fit to their combined 
data in Fig. [15] (red solid lines), showing that it 
underestimates the star formation rates we report 
here at both extremely high and low gas column 
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densities. 

At high column dens ities, the shallow portion of 
the bent power law of iHeiderman et al.l ( 20101 ) is 
driven by the inclusion of the I/MS'-HCN-derived 
data in the fit; their Spiteer-derived data seem rel- 
atively consistent with our analysis (examined in 
more detail below). That work includes a dis- 
cussion of the likelihood that far-IR-derived star 
formation rates may be underestimated for in- 
dividual star-forming regions by up to an order 
of magnitude. We find plenty of evidence for 
that to be the case in the nearby embedded clus- 
ters, where we can do direct comparisons between 
the IRAS far-IR-derived star formation rate and 
that derived from the YSO counts. For exam- 
ple, using the kno wn far-IR luminosity (26,000 L Q ; 
Ridge et al.ll2003 ) and the YSO count for the clus- 
ter core of the MonR2 m ain cloud core (132 YSOs, 
0.73 pc effective radius; iGutermuth et al.l |2009() . 
we obtain star formation rate surface densities of 
3 and 20 Mq yr _1 kpc~ 2 , respectively. A simi- 
lar discrepancy can be found in the Orion Neb- 
ula Clus ter. Its far-IR lum inosity is estimated at 
10 5 1/0 (jGezari et al.lll9 98). yet it s central parsec 



conta ins of order 10 3 YSOs (e.g. iFeigelson et al 



20051 ). yielding SFR surface density estimates of 
20 and 300 M Q yr" 1 kpc -2 , respectively. In sum- 
mary, the common conversion of far-IR flux to star 
formation rate appears to underestimate the true 
star-formation rate at parsec scales in nearby re- 
gions. 

We are not concerned with the apparent 
discrepancy with our result at low gas col- 
umn densities. The p rotostellar number counts 
(JHeiderman et al.ll2010h at these gas column den- 
sities are very low: a star surface density of 
10~ 17 Mq kpc~ 2 yr" 1 corresponds to one 0.5 Mq 
protostar per 0.5 Myr (the estimated lifetime of 
the evolutionary phase) in the area of an entire 
typical molecular cloud (~ 10 pc size scale is typ- 
ical of nearby clouds). We used the fitting tech- 
nique described in Section [3] to fit their Spitzer- 
derived data, excluding the upper limits. We find 
that the Spitzer-only correlation to be a power law 
of slope 3.7 ± 1.6 (dashed red line in Fig. IT51) . a 
steeper power law than our reported trend (though 
within la) and lacking in an apparent gas surface 
density threshold for star formation. Differences 
at this level may be explained in part by differ- 
ences in source classification method and uncer- 



tainty in the relative durations of the pertinent 
YSO evolutionary stages, but detailed consider- 
ation of such effects is beyond the scope of this 
paper. 

5. Summary 

We have presented an analysis of YSO and gas 
surface density at 0.3 to 2.0 pc scales for over 7000 
YSOs in eight nearby molecular clouds from 40 
square degrees of Spitzer- and 2MASS-surveyed 
sky. These plots show the relationship between 
the surface density of stars formed over the life- 
time of the molecular cloud vs. the surface density 
of the remaining gas in the clouds. A power law 
correlation between these two quantities is evident 
in all eight molecular clouds included in the study. 
In summary: 

• MonR2 and Ophiuchus have particularly 
strong and well-defined correlations, thus 
we fit them each with power laws (indexes 
of 2.67 and 1.87, respectively). 

• Some clouds exhibit large deviations from 
the power law form apparent in MonR2 and 
Ophiuchus; we argue that these deviations 
are due to gas dispersal by massive stars and 
non-coevality in each cloud. The presence of 
non-coeval regions in the clouds is demon- 
strated by variations in the ratio of the num- 
bers of Class II to Class I YSOs. 

• Even given non-coevality and gas dispersal 
in most clouds, we find that greater than 
half of all YSOs with excess infrared emis- 
sion in every cloud considered are found to 
have local £*/5] 2 in the range of 3 x 10~ 4 
and 5 x 10~ 3 pc 2 M" 1 . 

• Extracting the regions of each cloud found 
within a common evolutionary state ( 3 x 
10~ 4 < £*/£ 2 as < 5 x 10~ 3 ), the Class II 
to Class I ratio is shown to decline with gas 
column density. However, the magnitude 
of the decline varies by cloud, from near- 
uniform ratios with small declines in MonR2 
and North America Nebula to steep declines 
in SI 40 and Perseus. 

• The correlation can be reproduced by a sim- 
ple model where the star formation rate per 
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unit area varies with the square of the sur- 
face density of gas. This model accounts 
for the gas removed from the cloud by the 
accretion and ejection of gas mass by each 
star. We find that the accretion and ejec- 
tion of gas will cause deviation from the 
power law form, but will only result in a mild 
steepening of the correlation at high gas sur- 
face densities over the likely star formation 
timescales of a few Myr. In contrast, we find 
star formation laws where the star formation 
rate per area varies linearly or with the 1.4 
power of the gas surface density are not able 
to fit the observed correlation. 

• The observed dependence of the stellar sur- 
face density with the gas density squared is 
a prediction of the thermal fragmentation 
of a sheet-like isothermal layer. To achieve 
the observed range of local column densities, 
one must invoke some sort of nonuniformity 
in the local column density structure of the 
sheet, such as a "hub-filament" geometry as 
described and analytically characterized by 
Mverd (|2009al ). 



The correlation reported here is largely 
consistent with recent Spiizer-derived esti- 
mates of star formation rate and gas sur- 
face densitie s both averaged over clouds 
( Evans et al. 2003) and mea sured within 



them ( Heiderman et al.l [2010J). We have 



identified some systematic effects that may 
explain differe nces between our result s and 
the analysis of iHeiderman et al.l ( 20101 ) that 



probes comparable size scales. Regardless of 
these differences, all of the work considered 
consistently suggests that extragalactic star 
formation rate and gas surface density cor- 
relations significantly underestimate those 
seen in nearby clouds. 
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A. Expansion of the c2d Extinction Maps 

The c2d extinction maps are generated using extinction measurements from each source classified as 
having photospheric colors in their merged 2MASS and Spitzer catalogs. The measurement is a parameter 
returned as part of a blackbody fit to the 1-8/im spectral energy distributions of those sources, limited by 
what detections where achieved to reasonable S/N in that range. Unfortunately, that means that the c2d 
extinction maps are limited to the field of view at least covered by IRAC. Since the c2d survey areas were 
chosen to match elevated extinction in the surveyed clouds (Ay > 3), obtaining an unbiased local baseline 
is impossible from the maps directly. This baseline is needed to remove the foreground extinction to the 
reddening map. 

In order to obtain a reasonable baseline Ay offset for the clouds surveyed by c2d, we had to expand those 
maps such that they extend beyond the IRAC-surveyed region, which was largely constrained to Ay > 3 
areas in each cloud. Our process is as follows: 



Const ruct an Ay map of a large area around each cloud using the technique from iGutermuth et al 
( 20091 ) that was used for the non-c2d clouds included in this work. 



• Resample the large map onto the c2d map's grid. 

• Fit the Ay(c2d) vs. Ay(G09) correlation between the two maps using those pixels that have valid 
data in both maps, in the range where both maps should have strong agreement: 2 < Ay(c2d) < 10. 

• Convert the large map's Ay values to use the c2d map's effective reddening law and zero-point cali- 
bration using the fitted parameters and fill in the empty parts of the c2d map's grid. 

The fits are shown in Figs. [T6l H3 & UHl and the results are summarized in Table |U The figures reveal 
a consistent deviation and spreading in the data for Ay(c2d) > 12 where the Ay(c2d) values are generally 
higher than the fit. The declining background star density with increasing extinction forces the adaptive 
smoothing of the G09 method to smooth over poorly sampled high column density structure, diluting it in 
favor of higher statistical S/N. In addition, the inclusion of the IRAC-derived data in the c2d maps likely 
enables higher S/N extinction measurements in greater numbers in intermediate extinction regions. 

The fit results are summarized in Table |U The inclusion of IRAC data in the c2d extinction measurement 
process could lead to deviations in the effective reddening law; indeed, color excess ratios in various com- 
binations of IRAC a nd 2MASS bandpasses have been shown to vary slightly (e.g. Serpens and Ophiuchus; 
Flaherty et al.ll2007l) . That said, these effects are expected to be small, and that is born out by the < 15% 
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variation indicated by the near-unity slopes of the fits. The offset calibration deviations are easily explained 
by the use of local reference fields in the c2d maps to establish the typical colors of field stars near each 
cloud. In constrast, the G09 technique adopts an intrinsic H — K = 0.2 color for field stars and then utilizes 
constant baseline subtraction away from the cloud to account for local variations. With our cross calibration, 
we ensure a fair baseline treatment for the c2d clouds, regardless of the differences in their construction. 

Additionally, we can test how the extinction mapping technique affects the log £* vs. log T, gas fit to the 
Ophiuchus cloud; the power law index is 1.93 ±0.03 with reduced x 2 of 5.01 using the map generated by our 
technique, and 1.87 ± 0.03 with reduced x 2 of 4.78 using the c2d map. Thus the values differ by less than 
10% and are within 1.5 sigma of each other. 

This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project 
of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute 
of Technology, funded by the National Aeronautics and Space Administration and the National Science 
Foundation. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. 
This research has made use of the VizieR catalogue access tool, CDS, Strasbourg, France. This work is based 
in part on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion 
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Laboratory, California Institute of Technology under a contract 1407 with NASA. Support for the IRAC 
instrument was provided by NASA through contract 960541 issued by JPL. 
Facilities: Spitzer. 
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Table 1 
Data and Derived Properties by Cloud 



Name Mass a (Mg) 


Area (pc' 


! ) 


Resolution ( 


pc) Baseline Ay 


Mean A v " 


Max. A v a 


Nil 


Nj 


Dist. (pc' 


) Disl 


.. R 


of. 


YSO Ref. 


MonR2 


25800 


1170.0 






0.95 




-0.5 


1.5 


17.1 


815 


188 


830 




1 




9 


CcpOB3 


16100 


551.0 






0.42 




0.9 


1.9 


21.6 


1976 


196 


700 




2 




9 


S140 


7150 


200.0 






0.51 




0.8 


2.4 


21.3 


470 


56 


750 




3 




9 


Orion 


33200 


572.0 






0.75 




-1.0 


3.9 


24.1 


2861 


502 


414 




4 




10 


NANeb 


24400 


548.0 






0.29 




1.1 


3.0 


22.0 


745 


425 


600 




5 




5 


Oph 


3070 


43.0 






0.20 




0.2 


4.8 


38.0 


179 


71 


150 




6 




ll b 


Perseus 


4300 


72.8 






0.29 




0.7 


3.9 


31.5 


243 


111 


250 




7 




ll b 


Serpens 


2590 


45.0 






0.48 




3.7 


3.8 


18.7 


141 


54 


415 




8 




ll b 


References. — 


- (1) ICarpenterl 120001); 


(2) IMoscadelli et al.l 12009); 


(3)IHirota et al.l (20081 


il); (4)IMenten et al.l 120071); 


(5) IRebull et al.1 1201 


.1); 


(6) IWilking et al.l 


(2003); (7)IEnoch et alj 1200 


61); (8)IDzib et al.1 120101); (9) 


IGutermuth et all lin 


prep.); (10) Mcgcath ct al. ( 


in prep J); (11 


) lEvans et al 


.1 120091) 








a Integrated in 


ass and mean and maximum Av values are derived from the b; 

count discrepancies between our reported counts and Table 5 ( 
advisement, we have included counts derived from the sourct 


asclinc-subtra 


,cted extinction 
J 120091) have be 


maps for the areas within the Spitzer cov 

!cn confirmed by the authors as a mistake 
jOSJ), and their Table 5 source counts sur 


oragc. 


Minor source 


if Evans et al 


in their tabic. 


Under the their 


; catalogs of |Evans et al.| 12C 


amary will be 



corrected in a forthcoming paper on all of the Gould Belt clouds (M. Dunham, private communication). 



Table 2 
log £* vs. log Y> gas Fit Coefficients 



Indc 



Calibration Offset Power Law Index Reduced \ 



MonR2 


0.87 


CopOB3 


0.17 


S140 


0.52 


Orion 


0.60 


NANeb 


0.67 


Ophiuchus 


0.83 


Perseus 


0.59 


Serpens 


0.83 



4.42 


2.67 ±0.02 


6.22 


2.21 


1.77 ±0.01 


20.3 


1.84 


1.37 ±0.03 


7.81 


2.40 


1.80 ±0.01 


17.9 


3.54 


2.15 ±0.02 


12.2 


3.10 


1.87 ±0.03 


4.78 


-7.0 


3.8 ±0.1 


13.0 


3.05 


1.95 


a 



a Thc Serpens fit's \ minimization did not converge using the Flaherty ct al. { 2007) fitting 
technique; a simpler, unweighted bisector fit has been used to obtain the reported coefficients. 
Saturation in the gas column density appears to cause a dual valued signal in the high column 
density region of the trend. 
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Table 3 

E*/£g OS Median Values by YSO Class 

Name Class II Class I 

MonR2 0.86 x 1CT 3 0.68 x 10~ 3 

CepOB3 3.0 x 10~ 3 1.0 x 10~ 3 

S140 1.7 x 10~ 3 1.2 x 10~ 3 

Orion 1.7 x 10~ 3 1.0 x 10~ 3 

NANeb 0.64 x 10~ 3 0.81 x 10~ 3 

Ophiuchus 0.47 x 10~ 3 0.40 x 10~ 3 

Perseus 1.9 x 10~ 3 0.90 x 10~ 3 

Serpens 0.76 x 10~ 3 0.57 x 10~ 3 
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Table 4 
Extinction Comparison Fit Coefficients 

Name Constant Offset (Ay) Slope 

Oph 0.666 0.959 ± 0.001 

Perseus 1.127 1.126 ± 0.004 

Serpens 1.633 0.922 ± 0.005 



22 




Fig. 1. — Extinction map of the MonR2 cloud 
overlaid in red with the spatial distribution of 
Spitzer-identified YSOs. The inverted grayscale 
is a linear stretch from Ay = — 1 to 10 mag- 
nitudes. Contour overlays start at Ay — 3 mag 
and their interval is 2 mag. The IRAC coverage is 
marked by the green boundary. The projected po- 
sitions of the YSOs in MonR2 closely trace almost 
all of the areas of detectably elevated extinction. 
Denser clusters of YSOs are clearly apparent in 
the zones of highest extinction. 




Fig. 2. — Extinction map of the CepOB3 cloud 
overlaid in red with the spatial distribution of 
Spitzer-identified YSOs. The grayscale and con- 
tour properties are identical to those in Fig. [1] 
As in that figure, YSOs are predominantly pro- 
jected on the elevated extinction zones within the 
cloud, and clusters are found in the highest ex- 
tinction zones. However, unlike MonR2, the large 
CepOB3b young cluster in the northwest corner 
of the coverage is largely offset from significant 
extinction. Focussed examination of this region in 
particular suggests that the OB stars present have 
dispersed much of the local natal cloud material 
(JGetman et al.ll2009t lAllen et al.llin prep.1 ). 




Fig. 3. — Extinction map of the S140 cloud over- 
laid in red with the spatial distribution of Spitzer- 
idcntified YSOs. The grayscale and contour prop- 
erties are identical to those in Fig. [TJ 




45 35 

Right Ascension 



Fig. 4. — Extinction map of the Orion cloud over- 
laid in red with the spatial distribution of Spitzer- 
idcntified YSOs. The grayscale and contour prop- 
erties are identical to those in Fig. [TJ 
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Fig. 5. — Extinction map of the NANeb cloud 
overlaid in red with the spatial distribution of 
Spitzer-identified YSOs. The grayscale and con- 
tour properties are identical to those in Fig. [TJ 




Fig. 6. — Extinction map of the Ophiuchus cloud 
overlaid in red with the spatial distribution of 
Spitzer-identiiied YSOs. The grayscale and con- 
tour properties are identical to those in Fig. [TJ 




Fig. 7. — Extinction map of the Perseus cloud 
overlaid in red with the spatial distribution of 
5piteer-identified YSOs. The grayscale and con- 
tour properties are identical to those in Fig. [TJ 




Rr^nt Ascension 



Fig. 8. — Extinction map of the Serpens cloud 
overlaid in red with the spatial distribution of 
Spiteer-identified YSOs. The grayscale and con- 
tour properties are identical to those in Fig. [TJ 




£„ (Mo PC"') 

Fig. 9. — YSO surface density versus gas col- 
umn density for the eight molecular clouds con- 
sidered here. YSO surface density is computed 
from the local NN11 surface density of YSOs sam- 
pled at the position of each source. Gas col- 
umn density is derived from the nearest sample 
to each YSO's position in the extinction map. 
Black hashed overlays mark out the regions of 
Ay < 1 (ie. consistent with zero, given typical 
uncertainties) and where the angular surface den- 
sity of sources is less than 7 per square degree 
(consistent wit h residual contam i nant sources ac- 
cording to the iGutermuth et al.l (|2009l ) classifica- 
tion scheme). The green polygon marks the fit 
to, and estimated spread in, the MonR2 cloud's 
data, and the blue polygon is the same, but for 
the Ophiuchus cloud's data. 
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Fig. 10. — Same as in Fig. |9j but now the stel- 
lar surface densities are sampled uniformly by po- 
sition on the same uniform grid as the extinc- 
tion map instead of centered on the sources them- 
selves. This area-sampled representation better 
represents regions of high gas column density but 
low YSO surface density. 




Fig. 11. — Histograms of £*/£g as using the data 
in Fig. [§J separated by YSO evolutionary class. 
The Class II YSO histogram is in black, the Class I 
YSO histogram is in red, and each is normal- 
ized to their peak bin. The green histogram is 
the Class II histogram, but normalized to the 
Class I histogram's peak bin and then divided by 
30. This represents the upper limit to the ex- 
pected frequency that edge-on C lass II sources will 



be misclassified as Class I YSOs (Guter muth et al. 



2009). Vertical dashed lines at T,„/T, 2 gas = 3 x 10~ 4 
and 5 x 10 -3 mark the empirical boundaries be- 
tween different apparent evolutionary distinctions 
in the cloud: young, protostar-rich regions are on 
the left, most of the embedded cloud populations 
are found in the center and exposed, regions that 
have undergone gas dispersal are on the right. 
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Fig. 12. — Class II to Class I source count ratio 
versus gas column density for those objects with 
3 x 10~ 4 < £*/£g as < 5 x 10~ 3 . Among all clouds, 
there is a clear trend toward lower ratios at higher 
gas column densities. This could suggest that low 
gas column YSOs are relatively more evolved if 
they formed in situ (which could be suggestive of 
either an earlier onset of star formation in the cur- 
rently low column zones or that the YSOs in high 
column zones have longer Class I phase lifetimes). 
Alternatively, the low column YSOs may have mi- 
grated from high density regions or may be the re- 
sult of gas dispersal. This latter explanation may 
be particularly relevant for Cep OB3b, S140 and 
Perseus where there is clear evidence of gas dis- 
persal. 
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Fig. 13. — Results from the simple gas depletion 
model for an assumed star formation law of power 
law index 2 and an adopted star formation rate 
coefficient k = 9.5 x 10~ 4 Myr _1 pc 2 A/" 1 . Gray 
lines trace the star-forming evolution of parcels of 
gas of a given surface density through its evolution 
predicted by the model; as mass is converted from 
gas mass to stellar mass, the initial samples drift 
upward. For tracks that approach the regimes 
where t ~ to, gas is significantly depleted by the 
star formation process and points arc to the left 
as they rise. The black lines are isochrones defin- 
ing the locus of points for the coeval evolution of 
molecular gas spanning a range of initial gas den- 
sities. The ages for each isochrone are noted, in 
Myr. The colored polygons demarcate the area 
that the observations of the MonR2 (green, steeper 
rise) and Ophiuchus (blue, shallower rise) clouds 
occupy. The dot-dashed lines mark where t = t a 
and t = 0.1 io for the model parameters chosen. 
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Fig. 14. — Several other realizations of the 
model; the left column assumes an initial power 
law of slope 1 and k = 0.1 Myr -1 , and the 
right column assumes a slope of 1.4 and k = 
0.05 Myr" 1 pc 08 Mq 0A ; c = 0.5 in all realizations. 
These values of k were adopted to overlap with the 
observations, and do not affect the shape of the 
isochrones. The black isochrone lines are adopted 
for the same time steps as used in Fig. [T3] The 
top row has no implicit gas column density thresh- 
old for star formation to occur. The bottom row 
includes such a threshold at 100 M pc~ 2 ; be- 
low that column density threshold, no stars are 
formed. The MonR2 and Ophiuchus data extent 
polygons (green and blue, respectively) and t = to 
and t = 0.1 to dot-dashed lines are overplotted as 
in Fig. Q2J The slope of 1 models are clearly in- 
consistent with the fits, and the slope of 1.4 with 
the star formation threshold is also inconsistent. 
The slope of 1.4 model with no threshold is con- 
sistent with the Ophiuchus fit if some degree of 
evolution is applied, as the approximate effect is a 
steepening of the trend with time. 
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Fig. 15. — Our reproduction of the Star Formation 
Rate Surfa ce Density v s . Ga s Column Density 
figure from lEvans et al.l (|2009h . The filled black 
circles are mean measurements for entire molecu- 
lar clouds from the c2d survey as presented in that 
work, and the open square is their mean. Red dia- 
monds are similar cloud-integrated values derived 
from the data compiled in Table [TJ The black 
lines are s ome recent correlations reported in the 
literature feennicutt et~aflll998t IWu et all 120051: 
Bigiel et aUl2008l ) based largely on extra-galactic 
measures (500 pc or larger size scales). Dashed 
portions of the line are extrapolated relative to 
the reported range of each trend. The green and 
blue polygons demonstrating the correlations re- 
ported for MonR2 and Ophiuchus, respectively, 
have been converted to star formati on rate den- 
sity u nder identical assumptions to lEvans et al.l 
(J2009I ) and overlaid. The cloud integrated values 
are consistent with our measurements, but the lat- 
ter extend to higher column density because of the 
smal ler size scales con s idered . The red trends are 
from lHeiderman et al.l ( 20101) ; the solid line is their 
bent power law fit to a combined set of Spitzer- and 
JiMS'-derived data, and the dashed line is our fit 
to their Spitzer-derived data only. 



Fig. 16. — Ay comparison between the c2d- 
generated map of Ophiuchus and our own, resam- 
pled to their grid. The fit was performed on the 
gray points only, where we expect both maps to 
have strong agreement. 
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Fig. 17. — Ay comparison between the c2d- 
generated map of Perseus and our own, resampled 
to their grid. The fit was performed on the gray 
points only, where we expect both maps to have 
strong agreement. 
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Fig. 18. — Ay comparison between the c2d- 
generated map of Serpens and our own, resampled 
to their grid. The fit was performed on the gray 
points only, where we expect both maps to have 
strong agreement. 
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